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Abstract 



We anafyze the renormalon ambiguities that appear in factorization formu- 
las in QCD. Our analysis contains a simple argument that the ambiguities 
in the short-distance coefficients and operator matrix elements are artifacts 
of dimensional-regularization factorization schemes and are absent in cutoff 
schemes. We also present a method for computing the renormalon ambiguities 
in operator matrix elements and apply it to a computation of the ambigui- 
ties in the matrix elements that appear in the NRQCD factorization formulas 
for the annihilation decays of S-wave quarkonia. Our results, combined with 
those of Braaten and Chen for the short-distance coefficients provide an 
explicit demonstration that the ambiguities cancel in the physical decay rates. 
In addition, we analyze the renormalon ambiguities in the Gremm-Kapustin 
relation and in various definitions of the heavy-quark mass. 
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I. INTRODUCTION 



In Quantum Chromodynamics (QCD), it is often useful to describe physical processes 
involving more than one distance scale by making use of a factorization formalism. In such 
a formalism, a physical observable is written as a sum of products of short-distance coeffi- 
cients with long-distance operator matrix elements. The short-distance coefficients may be 
calculated as a perturbative series in the strong coupling constant a s , evaluated at the short- 
distance scale. The operator matrix elements contain all of the sensitivity of the physical 
observable to low-momentum (infrared) processes. Because of this infrared (IR) sensitivity, 
the operator matrix elements are generally not amenable to a perturbative calculation and 
are usually determined by comparison of physical quantities with experimental values or 
through nonperturbative methods, such as lattice QCD. Some well-known examples of this 
approach are the light-cone expansion in deep-inelastic scattering, the factorization expres- 
sions for hard hadron-hadron cross sections, Heavy-Quark Effective Theory (HQET), which 
is useful in the study of heavy-light mesons, and Nonrelativistic Quantum Chromodynamics 
(NRQCD), which is useful in the study of heavy quarkonium. 

The conventional wisdom is that the perturbation expansions for the short- distance co- 
efficients in QCD factorization formulas are, at best, asymptotic series §. One can see 
why this might be the case by examining the behavior of the perturbation series in the 
context of a simple model, which gives the leading large-iVy behavior of the theory. This 
model has seen extensive applications in recent years in such topics as the Operator Prod- 
uct Expansions (OPE) for e + e~ annihilation and r decay [§-0], Heavy Quark Effective 
Theory (HQET) [f?HT0||, and the NRQCD factorization formalism [l]. In this model, which 
we call the "bubble-chain model," one generates the "perturbation series" by inserting into 
the gluon propagator in a one-loop diagram all numbers of order-a s fermion-loop vacuum- 
polarization corrections. The resulting perturbation series gives the order-a s fermion-loop 
contribution to the running of the coupling constant. In this expression, one then replaces 
the fermion-loop contribution to the one-loop QCD beta function (3$ with the full one-loop 
QCD beta function. As we will describe later, in the bubble-chain model, the perturbation 
series diverges, with terms growing as the factorial of the order. 

One can attempt to resum the series by carrying out a Borel transformation. In the Borel 
plane, the factorial growth of the series corresponds to singularities, known as renormalons. 
If the renormalon singularities lie on the positive real axis, then they make the inverse 
Borel transform ill-defined. These ambiguities in the inverse Borel transform are potentially 
important because they have the same nominal size as corrections in the QCD factorization 
formulas that go as powers of the ratio of scales. 

In a well-defined theory, physical observables, in contrast with short-distance coefficients, 
are unambiguous. Of course, if one were to attempt to compute a physical observable 
entirely in perturbation theory, without making use of a factorization formalism, then one 
might encounter ambiguities, owing to the failure of the perturbation series to converge. 
However, such ambiguities are not intrinsic to the observable. They are artifacts of the 
method of computation and would be absent if one were to compute the observable in terms 
of nonperturbative expressions, such as lattice path integrals. 

When one computes a physical observable by making use of a factorization scheme that 
is defined in perturbation theory, such as one that is based on dimensional regularization, 
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then the short-distance coefficients and operator matrix elements can contain ambiguities 
that arise in consequence of the defining perturbation series. It is these ambiguities that 
arise from the factorization scheme that are the focus of this paper. 

If one were to attempt to compute an operator matrix element entirely in perturbation 
theory, then additional ambiguities, which have nothing to do with the factorization scheme, 
might arise. Such ambiguities are not intrinsic to the matrix element, since they would not 
be present if one defined the matrix element nonperturbatively and evaluated it in terms of 
nonperturbative expressions. As we shall describe later in this paper, one can isolate the 
ambiguities that arise from the factorization scheme by computing the difference between 
an unambiguous definition of the operator matrix element (for example, a lattice definition) 
and an ambiguous definition (for example, one based on dimensional regularization) . In the 
remainder of this paper, when we refer to ambiguities in operator matrix elements, we mean 
only those ambiguities that arise from the factorization scheme. 

Because physical observables are unambiguous, the renormalon ambiguities in the short- 
distance coefficients must be cancelled by corresponding ambiguities in the operator matrix 
elements, to the level of accuracy of the factorization procedure. Factorization procedures 
that are based entirely on an underlying effective field theory, such as those in the light-cone 
expansion, OPE, HQET, and NRQCD, are, in principal, valid to all orders in an expansion 
in the inverse of the large scale in the process. In contrast, factorization procedures for 
hadron-hadron-induced hard-scattering processes, for example, the Drell-Yan process, are 
known to fail for contributions that are subleading in the inverse of the large scale [jnj. In 



the discussions in the remainder of this paper, we have in mind the factorization procedures 
that are based entirely on underlying effective field theories. However, our results may also 
apply to other factorization procedures at the level of leading orders in the inverse of the 
large scale. 

If one determines operator matrix elements by comparing factorized expressions for physi- 
cal observables with results from experiment, then the cancellation of renormalon ambiguities 
holds by construction, for those observables. On the other hand, if one determines operator 
matrix elements through a calculation in the underlying effective field theory, as we do in 
this paper, then the cancellation of renormalon ambiguities is a nontrivial confirmation of 
the accuracy of the effective field theory in reproducing the low-momentum behavior of the 
full theory and of the consistency of the definitions of the short-distance coefficients and 
operator matrix elements. 

In this paper, we investigate the relations between the renormalon ambiguities in the 
short-distance coefficients and the renormalon ambiguities in the operator matrix elements. 
In our analysis, it is clear that the short-distance coefficients and the operator matrix ele- 
ments display no ambiguities in a factorization scheme in which ultraviolet (UV) divergences 
in the operator matrix elements are controlled with a regulator, such as lattice regulariza- 



tion, that has an explicit cutoff in the loop momenta [12,13]. In contrast, in a dimensional- 



regularization scheme, both the short- distance coefficients and the operator matrix elements 
contain renormalon ambiguities, but they cancel in physical observables JTO] . We give a pre- 



scription for calculating the renormalon ambiguities in the dimensionally-regulated operator 
matrix elements, and we apply this method to the computation of the renormalon ambigu- 
ities in the NRQCD matrix elements that appear in the decays of S-wave quarkonia in the 
leading and first subleading orders in the expansion in the heavy-quark velocity. Compar- 
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ing with a calculation of the ambiguities in the short- distance coefficients for the decays of 
S-wave quarkonia 0, we find that the renormalon ambiguities cancel in the physical decay 
rates, as expected. 

The remainder of this paper is organized as follows. In Sec. [H], we review the bubble- 
chain model and Borel transformation and discuss, in general, the origins of renormalons 
in short-distance coefficients. Here we present a concise argument to show that renormalon 
ambiguities in the short- distance coefficients arise from the low-momentum regions of loop 
integrals. Sec. |II| contains a discussion of the interplay between the factorization scheme and 
the renormalon ambiguities. In this section we clarify the fact that renormalon ambiguities 
in the short- distance coefficients and operator matrix elements are artifacts of dimensional- 
regularization factorization schemes and are absent in cutoff schemes. In Sec. |IV|, we describe 
our method for computing the renormalon ambiguities in the operator matrix elements. We 
use this method, in Sec. |V|, to calculate the renormalon ambiguities in the NRQCD matrix 
elements for decays of S-wave quarkonia. In Sec. |V11| , we examine the renormalon content 
of the various quantities that appear in the Gremm-Kapustin ]14| relation. We also discuss 
the ambiguities in various definitions of the heavy-quark mass. Finally, we summarize out 
results and discuss their implications in Sec [VII]] . 



II. RENORMALONS IN THE BUBBLE-CHAIN MODEL 

As we have mentioned, the bubble-chain model, gives the exact large-Ay behavior of 
QCD. In it, one generates the "perturbation series" by inserting into the gluon propagator 
in a one-loop diagram all numbers of order-a s fermion-loop vacuum-polarization corrections. 
That is, one replaces the factor a s (n 2 ) in the one- loop diagram with 



« s (/i 2 ) £ [-/W/x 2 ) Hl 2 e c /f)] n (la) 

n=0 

= ^1 (lb) 

l + /W/i 2 )ln(Z 2 e c //i 2 ) 1 J 

= ««(Z 2 ), (lc) 



where / is the Euclidean momentum of the gluon in the one-loop expression, C is a 
renormalization-scheme-dependent constant (C = —5/3 in the MS scheme), (3q is the one- 
loop contribution to the QCD beta function, and we identify the expression (HD as the one- 
loop running coupling constant a^\l 2 ). At this stage, (3$ contains only the fermion-loop con- 
tribution. However, in the bubble-chain model, one includes the effects of gluons on the run- 
ning of the coupling by promoting (3q to the complete QCD expression /5 = (33/2— Nf) /(6ir). 
Now, suppose that a part of the original one-loop integral has the IR behavior 

A dir. (2) 

Here we have introduced a cutoff A on the magnitude of I, in order to focus on the low- 
momentum region. We assume that the integral is IR finite (m > —1), since any infrared 
divergences are absorbed into the operator matrix elements in the factorization procedure. 
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If one inserts the expression (|la|) into the one-loop integral (0) and integrates term by term, 
then it is easy to see, using 



I 



" dl l m ln n (/ 2 //i 2 ) = 2> m+1 (-l) n n!/(m + l) n+1 , (3) 



that terms in the resulting series grow as n\. Hence, the series is only asymptotic. 

One way to associate a well-defined function with an asymptotic series is through the 
Borel transform. If a function f(a s ) has a power-series expansion 

oo 

fM = a "«" > ( 4 ) 

n=0 

then its Borel transform is defined by 

OO 1 

f{t) = a S(t) + £ — — aX- X • (5) 

In this paper, we indicate the Borel transform of a quantity / either as / or B[f]. Clearly, 
the series (|5]) for the Borel transform has better convergence properties than the original 
power series (f|). The function f(a s ) that generates the series (H) can be recovered from the 
Borel transform f(t) through the inverse Borel transform 

f(a s ) = / dt e- 1 ^ f(t) , (6) 



provided that the integral (§) is sufficiently convergent and that the Borel transform (||) is 
well-defined along the positive real axis. [] 
The Borel transform of the series fllaD is 



3f '(«) = 

where u = Pot. If we insert the Borel transform (0) into the one- loop integral (g), we obtain 

„A \m+l-2u 

(^ 2 /e c r dll m - 2u = (^/e c r — - , (8) 
jo m + 1 — 2u 

where A is an ultraviolet cutoff, and, in the spirit of dimensional regularization, we have 
defined the integral for Re u > (m + l)/2 by analytic continuation. The integral has a pole 
in the complex u plane at u = (m+ l)/2. This singularity appears even though the original 
one-loop integral is perfectly well-behaved. A pole in u, such as this one, that is associated 
with the running of a s is called a renormalon 0. 

In order to understand the origins of this renormalon singularity more fully, let us re- 
examine the integral (H). Applying an IR cutoff A, we have 



1 If the power series (Q) is only asymptotic, then there is not a unique function f(a s ) that generates 



it. The inverse Borel transform picks out a function with particular analyticity properties [15] 
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«A Am+l-2u \m+l— 2u 

{f/e c Y dir-^ = ^/e c T- -— . (9) 

Jx m + 1 — 2u 

If A and A are finite positive numbers, then the integral (^) is well defined, even at the point 
u = (m + l)/2. Furthermore, the inverse Borel transform is convergent, provided that 

a s (jj 2 )(3o\n[/j 2 / (X 2 e c )] < 1, (10) 

which is identical to the condition that the integration over I never passes through the 
Landau pole in ([Td]) . Hence, we conclude that the renormalon singularity in the u plane 



arises from the region of integration near / = 0, which we have excluded by introducing the 
IR cutoff A. The renormalon singularity is a signal that, because of the growth of the running 
coupling, the perturbation series breaks down near zero loop momentum. A renormalon that 
arises from the region near zero loop momentum is called an IR renormalon. 

The renormalon in (|8|) is on the positive real axis. Hence, the inverse Borel transform is 
not unique. One can obtain various values, depending on the prescription used to deform the 
integration contour near the pole. For example, one can deform the integration contour in 
(^) into the complex plane so that it runs above the pole or below the pole, or one can take 
a linear combination of these two contours, such as the principal value. These prescriptions 
all differ by amounts that are proportional to the residue of the integrand of @ at the pole. 
If the pole in f(t) is at the point t* = u*/j3o and has residue R*, the ambiguity in f(a s ) has 
the form 

Af(a s ) = K{27TR*)e- tt/as ^ 2) , (11) 

where K is a constant of order unity. In a s (fi 2 ), the scale \i is typically chosen to be the 
largest scale in the physical process. For large /z 

A) ln(/x 2 /A4 CD ) 

Inserting this expression into (|ll|); we see that the renormalon ambiguity is given approxi- 
mately by 

Af « K(2nR*) (^jjrj • (13) 

Thus, the ambiguity is suppressed as a power of fi 2 . However, such an ambiguity can be 
of practical importance, since it is of the same nominal size as power corrections in the 
factorization formulas. 

We have seen that IR renormalons can arise from the region of integration near zero 
loop momentum. In an analogous fashion, UV renormalons can arise from the region of 
integration near infinite loop momentum. The loop integration extends to infinity in unreg- 
ulated, UV-finite integrals and in UV-divergent integrals in dimensional regularization. In a 
renormalizable theory, the loop integration in a short-distance coefficient is, at worst, loga- 
rithmically divergent in the UV region when u = 0. (The integrand that we have displayed in 
the expression (|8|) is an approximate form that is valid only in the IR region.) In the case of 
a logarithmic UV divergence, the pole at u = is identical to the dimensional-regularization 
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pole in e = (4 — D)/2, where D is the number of space-time dimensions, and it is removed 
in the standard minimal-subtraction renormalization procedure. The remaining pieces of 
the integrand, which have a convergent power count in the large-momentum region when 
u = 0, may contribute additional renormalons. However, by virtue of their convergent power 
count, these pieces yield renormalons only on the negative real u axis. Renormalons on the 
negative real u axis do not introduce ambiguities into the inverse transform in Eq. (|6|). 

Another potential source of renormalons in the short- distance coefficients is the UV 
renormalization procedure that is applied to full QCD. If the renormalization counterterms 
involve loop integrations down to zero momentum, then the counterterms can contain renor- 
malons. In this paper, we assume that a minimal-subtraction procedure (MS) or modified 
minimal-subtraction procedure (MS) is used, so that the renormalization counterterms are 
poles in e or constants and contain no IR renormalons. 

III. FACTORIZATION-SCHEME DEPENDENCE OF RENORMALONS 

In a factorization formalism, the domain of integration of a loop integral in full QCD 
is partitioned between the operator matrix elements and the short- distance coefficients. 
This partitioning allocates the low-momentum part of the loop integral, including the IR 
divergences, to the operator matrix elements and allocates the remaining high-momentum 
part of the loop integral to the short-distance coefficients, which contain no IR divergences. 

The partitioning arises in the matching of the effective theory to full QCD. For example, 
in the on-shell matching procedure, on-shell quark and gluon amplitudes in the full theory 
are equated to the same amplitudes in the effective theory. Each amplitude in the effective 
theory can be written as a sum of products of short-distance coefficients with matrix ele- 
ments in the on-shell quark and gluon states. Low-momentum parts of the amplitudes in 
the full theory, including all of the IR divergences, correspond to matrix elements in the ef- 
fective theory. These matrix-element contributions can be factored from the full amplitudes. 
Then, the remaining IR-finite parts of the full amplitudes correspond to the coefficient func- 
tions. In the perturbative implementation of the matching procedure, which is familiar, 
for example, from treatments of the QCD-improved parton model, the matrix-element con- 
tributions are subtracted from the full-QCD amplitudes order-by-order in a s . Then, the 
remainders from this subtraction procedure are identified with perturbative contributions 
to the short-distance coefficients in the effective theory. 

The amount of the full-QCD amplitude that resides in the matrix elements and, hence, 
in the short- distance coefficients, is controlled by a UV regulator that is imposed on the 
operator matrix elements.^ The UV cutoff plays the role of a factorization scale, and the 



2 Throughout this paper, we speak of the short-distance coefficients as being dimensionally regu- 
lated, lattice regulated, etc. What we mean is that these are the short-distance coefficients that 
correspond, in the factorization formalism, to operator matrix elements that are dimensionally 
regulated, lattice-regulated, etc. Of course, one can apply an IR regulator to the short-distance 
coefficients as an intermediate step in calculations. However, dependences on such IR regula- 
tors cancel in the matching conditions, and the IR behavior of the short-distance coefficients is 



7 



choice of UV regulator is known as the "factorization scheme." The consistency of the 
underlying effective field theory guarantees that physical quantities are independent of the 
factorization scheme, provided that one works to sufficient accuracy in the effective field 
theory. 

One class of UV regulator consists of those regulators in which the cutoff is manifest 
and alters the form of the integrand at large momentum to make it more convergent. We 
call such regulators "cutoff" regulators. Pauli-Villars regulators, lattice regulators, and 
explicit bounds on the magnitude of the Euclidean momentum are all examples of cutoff 
regulators. For such regulators, power UV divergences manifest themselves as explicit power 
dependences of the matrix elements on the UV cutoff. For our purposes, the crucial property 
of cutoff regulators is that they do not alter the behavior of the integral below the cutofffj 

In contrast, in dimensional regularization in D dimensions, the cutoff enters only implic- 
itly through an overall factor that goes to unity when 2e = 4 — D vanishes. After one discards 
poles in e, one finds that logarithmic divergences are cut off by the dimensional-regularization 
scale.0 However, terms with a power-divergent power count that are homogeneous in the 
integration momentum vanish.^ Because of this property, dimensional regularization poten- 
tially alters the behavior of the integral below the cutoff. 

In a cutoff factorization scheme, the loop-integration for the IR-finite parts of the short- 
distance coefficients does not extend to zero momentum. The part of the integration below 
the cutoff resides in the operator matrix elements, to the extent that the effective theory 
accurately reproduces the low-momentum behavior of the full theory. Hence, in a cutoff 
factorization scheme, the short-distance coefficients are free of a given IR renormalon am- 
biguity, provided that the effective action retains enough powers of the small scale in the 
calculation to reproduce the low-momentum behavior that is characteristic of such a renor- 
malon. The absence of renormalons in a cutoff factorization scheme has been discussed in a 



number of publications. Some of the earliest discussions are in Refs. ||12| , [13|1 . 

In a dimensional-regularization factorization scheme, it is conventional to integrate the 
IR-finite parts of the short- distance coefficients down to zero momentum. Consequently, 
the short- distance coefficients contain renormalon ambiguities. The IR-finite parts of the 



controlled by the UV regulator of the operator matrix element. 

3 Lattice regulators and Pauli-Villars regulators do alter the behavior of the integral below the 
cutoff, but only through terms that scale as inverse powers of the cutoff. In an effective field 
theory, one compensates for this effect, restoring the behavior of the theory below the cutoff to 
the required level of accuracy, by including in the effective action corresponding terms that scale 
as inverse powers of the cutoff. 

4 When we speak of dimensionally-regulated matrix elements in this paper, we assume that such 
poles, and possibly some associated constants, have been discarded. That is, by dimensionally- 
regulated matrix elements, we really mean MS or MS matrix elements. 

5 We use the term "homogenous" to denote expressions that are a pure power or a pure power 
times logarithms. 
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short- distance coefficients correspond to the power-UV-divergent parts of the operator ma- 
trix elements. Thus, because of the constraints imposed by the matching conditions, the 
integration of the IR-fmite parts of the short-distance coefficients down to zero momentum 
demands, for consistency, that the UV-power-divergent parts of the matrix elements be set 
to zero. As we will explain in Sec. |1V B , there is a prescription for regulating the operator 
matrix elements dimensionally in which such power UV divergences are set to zero. This 
prescription implies that the dimensionally-regulated matrix elements do not completely 
reproduce the behavior of the full theory below the factorization scale. 

Matrix elements are completely determined, in principle, in terms of short-distance coef- 
ficients and physical observables. Therefore, the absence of ambiguities in the cutoff short- 
distance coefficients implies that the cutoff matrix elements are also ambiguity-free. This is 
consistent with the fact that lattice-regulated matrix elements have an unambiguous defini- 
tion in terms of path integrals in the effective theory. In a cutoff scheme, the low-momentum 
content of the theory resides completely in the operator matrix elements. 

In a dimensional-regularization scheme, the short- distance coefficients contain ambigui- 
ties. Therefore, the dimensionally-regulated operator matrix elements must contain ambigu- 
ities that cancel those in the short- distance coefficients. The appearance of renormalon am- 
biguities in the short-distance coefficients and operator matrix elements in the dimensional- 
regularization scheme has no physical significance: It is a factorization-scheme artifact. 
We conclude that the factorial growth in the perturbation series for the short-distances 
coefficients is cancelled by the factorial growth in the perturbation series that relate cutoff- 
regulated matrix elements to dimensionally-regulated matrix elements, provided that one 
works to the same order in a s in both quantities, and provided that one works to sufficient 
accuracy in the underlying effective field theory [10]. 

Let us explain more fully what we mean by sufficient accuracy in the underlying effective 
field theory. Generally, an operator matrix element scales, in its leading behavior, as a 
power of the small scale in the two-scale problem. For example, in the NRQCD factorization 
formalism for heavy-quarkonium decay and production [16j , matrix elements scale as powers 
of p = mv, where p and v are the typical heavy-quark-antiquark relative momentum and 
velocity, and m is the heavy-quark mass. To achieve the cancellation of a renormalon at 
given value of u, it is crucial to retain matrix elements of sufficiently high order in the small 
scale to reproduce the IR power behavior that gives rise to that renormalon. Furthermore, if 
one is computing a matrix-element renormalon ambiguity in the context of an effective field 
theory, it is necessary to retain terms in the effective action of sufficiently high accuracy 
in the small scale to reproduce the corresponding IR power behavior. The action of an 
effective theory may depend on both the cutoff and the large scale. For example, in HQET 
and NRQCD, the action depends on m, as well as the cutoff. Hence, the small scale may 
appear in the ratio of the small scale to the large scale or in the ratio of the small scale to 
the cutoff. In such cases, it is necessary to retain terms of sufficient accuracy in both the 
expansion in inverse powers of the large scale and the expansion in inverse powers of the 
cutoff. 
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IV. METHOD FOR COMPUTING RENORMALON AMBIGUITIES IN MATRIX 

ELEMENTS 



A. Short-distance expression for the ambiguity 

To identify the renormalon ambiguities in a dimensionally regulated matrix element, we 
make use of the fact that corresponding cutoff- regulated matrix element is free of ambiguities. 
Then, instead of computing the dimensionally-regulated matrix element itself, we compute 
only the effects of a change of regularization scheme from cutoff to dimensional. 

The effects of a change of regularization are contained in finite renormalizations of the 
operators in the underlying effective theory. In general, in order to work out these finite 
renormalizations, it is necessary to identify all of the renormalization counterterms in the 
effective theory and to compute their finite coefficients. However, at the one-loop level, 
this amounts merely to computing the difference between the dimensionally-regulated and 
cutoff amplitudes. (In both the dimensionally-regulated and cutoff amplitudes, we use the 
dimensionally-regulated expression for the vacuum-polarization insertions in the gluon prop- 
agator.) 

Because the renormalization parts of an operator matrix element are short- distance quan- 
tities, we can analyze them perturbatively.0 Furthermore, the renormalization parts are 
independent of the external states. Therefore, we are free to make a choice of the external 
states that is convenient for a perturbative analysis, namely, states consisting of on-shell 
elementary quanta (quarks, gluons, etc.). 

We note that, for purposes of computing the ambiguities in the dimensionally-regulated 
matrix elements, any reference cutoff regulator will do. At the one-loop level, a particu- 
larly convenient choice of cutoff regulator is a simple cutoff on the magnitude of the three- 
momentum of the gluon in the loop. While such a cutoff is not, in general, consistent with 
gauge invariance, it is compatible with the QED-like nature of the gauge invariance at one 
loop. With this choice, we compute 

/•dim p r\ r rdim 

(/C) dim - (/C) cutoff = J q dH I K - J dl J o d?l I K = J dl J x dH I K , (14) 

where /C is an operator, 1% is the integrand corresponding to its matrix element between 
on-shell states in one-loop perturbation theory, and "dim" and "cutoff" denote dimensional 
and cutoff ultraviolet regularization, respectively. We write "dim" as the upper limit of 
integration to indicate dimensional regularization of UV divergences. A is a cutoff on the 
magnitude of the 3-momentum / and corresponds roughly to the IR cutoff A discussed in 
Sec. ||. We assume that A is sufficiently large that the condition (|l(J) is satisfied, so that 



the integration on the right side of Eq. ( [HD does not pass through the Landau pole. 

The expression (|T4]) , being a renormalization-scheme dependence, is IR finite, as is ap- 
parent from the form on the right side of the equation. However, it may happen that 



6 In fact, the dimensionally-regulated matrix elements have only a perturbative definition, either 
in terms of physical observables or in terms of matrix elements in some other scheme, such as 
lattice regularization, that does have a nonperturbative definition. 
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the separate terms in the left and middle expressions in Eq. ([14]) require an IR regulator. 
Then, one must, for consistency with the complete IR- finite expression, impose the same IR 
regulator on both terms. 

The expression on the right side of Eq. ([14]) is the basis for our calculation of the renor- 
malon ambiguities. One could equally well use the middle expression in Eq. (|T4|), taking 
advantage of the fact, explained in Sec. [IV B| , that the dimensionally- regulated integral van- 
ishes for IR-finite terms in the integrand, so that one need compute only the cutoff integral.]] 
However, we wish to emphasize the insensitivity of the renormalon ambiguities to the low- 
energy physics by working with the right side of Eq. (|nj), which contains no low- momentum 
contribution. With this computational procedure, the dimensionally-regulated matrix ele- 
ments contain only UV renormalons. 



B. Dimensional-regularization prescription 

Let us now explain more precisely what we mean by dimensional regularization of oper- 
ator matrix elements. It is well-known that dimensionally-regulated expressions can some- 
times be ambiguous when the integrand involves a limiting procedure. For example, the 
result for a dimensionally-regulated integral can be changed by carrying out formal manip- 
ulations involving power series in the integration variable. Consequently, it is necessary to 
specify how the integrand is to be arranged before the dimensional regularization is imposed. 

In computing the short- distance coefficients in the dimensional-regularization scheme, it 
is conventional to integrate IR-finite terms down to zero loop momentum. As we have dis- 
cussed in Sec. [TTT], such an approach corresponds to setting to zero the UV power divergences 
in the operator matrix elements. Therefore, for consistency with conventional calculations 
of the short- distance coefficients, we wish to follow a procedure for the dimensional regular- 
ization of the operator matrix elements in which UV power divergences are set to zero. Such 
a procedure is usually employed in regulating, for example, the matrix elements for parton 
distributions [p"7jl . 

One applies this dimensional regularization procedure to QCD corrections to the matrix 
elements of operators between on-shell quark and gluon states. Such expressions appear, for 



example, in Eq. (|Tj), and also in the matching conditions between full QCD and the effective 
field theory that fix the short-distance coefficients. One first decomposes the QCD correc- 
tions to the operator matrix element into a linear combination of the tree-level operators in 
the effective theory. (In the case of NRQCD, this amounts to Taylor expanding the integrand 
with respect to the external momenta and taking appropriate linear combinations to form 
quantities of definite orbital angular momentum, spin, and color.) All of the dependence on 
the external momenta now resides in the tree-level operators. The coefficients of the tree- 
level operators contain the integration over the loop momentum I and are now independent 



7 It is interesting to note that the cutoff integral in the middle expression in Eq. (|i~4| ) contains 
IR renormalons, whereas the right side of Eq. (|l4|) contains UV renormalons. Thus, we see that 
UV renormalons can be replaced by IR renormalons through the manipulation of quantities that 
vanish in dimensional regularization. 
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of the external momenta. In some effective theories, such as NRQCD, the coefficients of 
the bare operators may not be homogeneous in I. In that case, one further expands the 
integrand in a power series in / divided by the large scale (m in the case of NRQCD). Then, 
each term in the series is homogeneous in I. Now, any terms that are IR finite (that is, UV 
power divergent) require no IR regulator. Hence, such terms are scale invariant, and their 
integrals over all I vanish in dimensional regularization. We note that this formal procedure 
is equivalent to absorbing the UV-power divergences into a renormalization of the operator 
matrix elements. 

In light of these general considerations, we see that the application of Eq. (|T5p requires 
some care. To obtain a result that is consistent with conventional calculations of the short- 
distance coefficients in the dimensional-regularization scheme, it is essential to use a proce- 
dure for computing the dimensionally-regulated matrix elements, such as the one that we 
have specified, in which the UV power divergences vanish. 



V. NRQCD FACTORIZATION FORMALISM 

In order to establish our notation and conventions, we now give a brief review of the 
NRQCD factorization formalism for heavy- quarkonium decays. 

In the nonrelativistic limit, in which the heavy-quark (or antiquark) velocity v is much 
less than unity, full QCD is described accurately by the effective theory NRQCD. In particu- 
lar, NRQCD is convenient for describing heavy-quarkonium bound states. In such systems, 
there are three important energy scales: the heavy-quark mass ~ m, the heavy-quark 3- 
momentum ~ mv, and the heavy-quark energy ~ mv 2 . The NRQCD effective theory is 
constructed from full QCD by integrating out the effects from energies of order m or larger. 
Since NRQCD describes only the effects at energy scales less than m, there is no heavy- 
quark- antiquark pair creation in the effective theory. That implies, to first approximation, 
that the heavy-quark and heavy-antiquark fields are decoupled. Hence, it is convenient to ex- 
press the effective action in terms of Pauli two-component spinors — one for the heavy-quark 
field and one for the heavy-antiquark field. 

The NRQCD effective Lagrangian can be used to reproduce the full theory to any desired 
accuracy in v. We can write it, up to terms of relative accuracy v 3 , as |18| , |19| 



where 



£nrqcd = Co + 5£ , (15) 



C = V 1 ( iD + ) iP + X ] (iD - ^ X (16) 




is the leading-order NRQCD effective Lagrangian, and 



^ = £ I [^(D 2 )^-X t (D 2 ) 2 x] 
c 2 



c 3 



[^(D -gE-gE- D)V + X f (D ■ gE - gE ■ D) X ] 



+ ^r[^ f (iD x gE-gE x z'D) ■ (Tip + xHi~D x gE - gE x zD) • ax] 

8m 2 

+ -?L[^(gB-<T)lP-x\gB-CT) X ] (17) 
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contains the relative order- 1> 3 corrections. Here, = <9 M + igA^, if> is the two-component 
Pauli field that destroys a heavy quark, x is the two-component Pauli field that creates a 
heavy antiquark, and c\, c 2 , c 3 , and c 4 are dimensionless short- distance coefficients. They 
are determined by matching amplitudes in the effective theory with those in full QCD. At 
tree level, c\ = c 2 = C3 = C4 = 1. There is a spin symmetry in the lowest order Lagrangian 
(fL6|), since it is independent of the spin of the heavy quark and heavy antiquark. This 
symmetry is violated by the terms in the higher-order corrections ( |TTD that contain er. 

In heavy-quarkonium annihilation decays, several distance scales are involved. The an- 
nihilation of the heavy quark and antiquark occurs at the short- distance scale 1/m, whereas 
the dynamics of quark- antiquark binding involves, principally, the long-distance scales 1/mv 
and 1/mv 2 . The NRQCD factorization formalism fl6j separates the physical effects that oc- 



cur at the scale 1/m from those at the longer- distance scales. Since effects from energies of 
order m are integrated out in obtaining the effective theory, the details of the annihilation are 
not described in NRQCD. However, the amplitude for a heavy quark and antiquark first to 
annihilate and then to be recreated is taken into account through effective local four-fermion 
interactions. The coefficient of a four-fermion operator contains the information about the 
effects at the scale m; the matrix element of the four-fermion operator in a quarkonium state 
contains the information about the long-distance effects. 

The imaginary parts of the four-fermion interactions are related, through the optical 
theorem, to the total decay rate. Hence, the heavy-quarkonium annihilation decay rate can 



be written in the factored form [16 



T{H) = E C mn(H\Omn\H) , (18) 

where C mn is a short- distance coefficient that is proportional to the imaginary part of the 
coefficient, in the effective action, of the four-fermion operator, and Mn is the mass of 
the state H. The matrix elements (H\O mn \H) are expectation values of local 4-fermion 
operators in the quarkonium state H. These local operators have the general structure 

O mn = 1p ] K, m X X'fcn.lp , (19) 

where fC n and K m are direct products of a color matrix (1 or T a ), a spin matrix (1 or a 1 ), 
and a polynomial in the gauge- covariant derivative D and in the field strengths E and B. 

If K m and K n are color-singlet operators, one can use the vacuum-saturation approxi- 
mation [|1J to simplify the matrix elements. This approximation is obtained by inserting 
a complete set of intermediate states between the quark-antiquark bilinears and retaining 
only the vacuum hadronic state. The result is 

(H\O mn \H) « {H\^K mX |0)(0| x^n^\H) . (20) 

Hence, a matrix element can be factored into a product of the matrix elements (H \^K mX |0) 
and (0| x*fcnip\H). In the case of hadronic decays, this approximation is valid up to cor- 
rections of relative order v A |[L6|| . For electromagnetic decays, the insertion of a vacuum 
projection operator is exact, since there are no hadrons in the final state. 

Generally, the NRQCD matrix elements are nonperturbative in nature. However, it is 
important to realize that they are fully determined, in principal, by the NRQCD effective 
theory. In fact, they can be computed by making use of the lattice formulation of NRQCD 
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VI. COMPUTATION OF THE U = 1/2 RENORMALON IN NRQCD MATRIX 

ELEMENTS 

In this section, we apply the method described in Sec. [TV] to compute the u = 1/2 renor- 
malons, which yield the ambiguities that are leading in Aqcd / for the hadron-to- vacuum 
matrix elements (0|xVV|^/V , )> (^Ix^lVc) , (Olx^D 2- ^! J/ip)/m 2 , (0|x t D 2 ^|^ c )/^ 2 - These 
are the matrix elements that arise in leading and next-to-leading order in v 2 in the decays 
of the S-wave charmonium states. Following Ref. M, we employ a shorthand notation for 
these matrix elements: 



(KU = e-(0|xW>|JMe)>, (21a) 

(£)„ = (0|xVl*fc>, (21b) 

{K, D ^ = (l/m 2 )e-(0\x^B 2 ^\J/^(e)), (21c) 

{K D *) V = (1/m 2 ) (Oix^Vkc) • (21d) 



As we have discussed in Section |V|, the renormalon ambiguity of the NRQCD matrix 
elements [Eq. (|TJ])] is insensitive to the low-energy behavior of the theory. It is, therefore, 
independent of the external states. Thus, we can evaluate the operator matrix elements in 
free cc states, using perturbative NRQCD (pNRQCD). That is, we use pNRQCD to evaluate 
the matrix elements 



(JC) V = e-(0|xW|cc(e)>, (22a) 

(K) P ee (0|xV|cc), (22b) 

(K D z) v ee (l/m 2 )e- (0| X VD^|cc(e)), (22c) 

(K D *) P ee (1/m 2 ) (0| X + DV|cc) • (22d) 



In each state, we take the quark and antiquark to be at rest and on their mass shells. The 
latter choice makes the calculation manifestly gauge invariant. 

A. Feynman rules 

In the pNRQCD calculation, we take for the unperturbed Lagrangian the noninteracting 
part of leading-order- m-v NRQCD Lagrangian (|16[) . The interactions in Eq. (|16D , as well 
as the higher-order-in-f corrections in Eq. (|TTD are treated as perturbations. The Feyn- 
man rules may be derived by standard methods. The rules for the quark and antiquark 
propagators and the "Abelian" parts of the quark-gluon interaction vertices are shown in 
Table Q. By "Abelian" we mean the terms that contain no color-matrix commutators. For 
our computation, which is at the one-loop level, the non-Abelian terms do not contribute. 

We also need the rules for the gluon propagators with all numbers of fermion-loop vacuum 
polarization insertions. Each renormalized fermion-loop vacuum-polarization yields a factor 

HV" + m v ) , (23) 

where I is the gluon momentum. (As we have already mentioned, (3q contains only the 
fermion-loop contribution at this stage, but eventually we promote /3q to the complete QCD 



if3 a s 



In 



C 
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expression (3q = (33/2 — iV/)/(67r).) For each vacuum-polarization insertion in the gluon 
propagator, there is a factor (^) and a factor from an additional free-gluon propagator. 
The sum over all numbers of such vacuum-polarization insertions is a geometric series, 
which is easily computed for the gluon propagator in a given gauge. In a covariant gauge, 
the propagator with all numbers of vacuum-polarization insertions is 



oo 
n=0 



hi 



-I 2 -It 
M 2 



c 



l 2 + ie 



I 2 + ie ' 



(24a) 



where we have suppressed color indices. Similarly, in the Coulomb gauge, the Coulomb-gluon 
propagator is 



oo 

iD c = £ Po a? 

n=0 



In 



IE 



C 



,2 > 



(24b) 



and the transverse-gluon propagator is 



n=0 



hi 



-/ 2 



C 



i(g - rv/i 2 ) 

l 2 + ie 



(24c) 



In each instance, the a° term is the free-gluon propagator. Here, and throughout this 
paper, we use Greek letters for Minkowski-space indices and Roman letters for 3-space 
indices. We do not distinguish between upper and lower 3-space indices; both correspond 
to the upper Minkowski-space index in covariant quantities and the lower Minkowski-space 
in contravariant quantities. 

It is now straightforward to obtain the Borel transforms of a s times the propagators in 
Eqs. flU). They are 



iB[a.DZ] 



i(- g ^ + m v /i 2 ) , m v j\ 2 



(_/2 _ , e )l +l 



-I 2 - ie 



(25a) 



iB[a s D 



c 



:,C 



l Z (-P 



(25b) 



iB[a a D%\ 



- klj/l 2 ) 
e c (-l 2 -ie) l + u 



(25c) 



B. The computation 



Now we apply the method of Sec. |TV] to calculate the u — 1/2 ultraviolet renormalons 
associated with the matrix elements in Eq. (p2|). For each perturbative correction to these 
matrix elements, we compute the right side of Eq. flT4T), using the convention of Sec. [IV B| for 
the form of the integrand in dimensional regularization. We consider only the ambiguities 
that are proportional to the lowest-order matrix elements (JC)y and {K)p- To obtain the 
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mixing into these matrix elements, we set the external momentum equal to zero in pertur- 
bative corrections. We find it most convenient to carry out the calculation in the Coulomb 
gauge. 

First we identify the u = 1/2 ultraviolet renormalons associated with the matrix elements 
(k>D 2 )v and {Kd 2 )p- Consider the vertex corrections to these matrix elements that arise from 
a Coulomb-gluon exchange with two Aq vertices. These corrections contain pieces propor- 
tional to the lower-order matrix elements (JC)y and (JC)p, respectively, which are obtained 
by setting the external momentum equal to zero. Since the A vertices respect the NRQCD 
spin symmetry, the vector and pseudoscalar renormalization constants of proportionality are 
equal. They are given by 

cy cy , r iVV ^ m dH 1 f 1 

dZ V}D 2 = dZ P}D 2 = A-kiCf / -j--r7 - 2// : . j ~i 



(2tt) 4 l - I 2 /(2m) + it m 2 -l - I 2 /(2m) + it 
1 

I 2 (-P - ieY 



X 727 ,! ...w . (26) 



where J dim d 4 l is shorthand for J dlo J^ im d 3 l, and Cf = 4/3. To evaluate this integral, we 
first integrate over Iq, using the contour-integral method. It is easy to see that the cut 
contribution has a power count such that it cannot contribute a pole at u = 1/2. The 
contribution from the quark (or antiquark) pole yields 

A7 S7 167F />T dH 1 (97\ 

5Z v , D2 = 5Z P , D2 ~-—^-j __ f+2u _ _ [2/ ^ u . (27) 

Now we expand the factor [1 — I 2 / (4m 2 )] -1 in a power series in Z 2 /m 2 . Only the first term 
in the series has the correct power count to yield a pole at u = 1/2. Integrating over I, we 
find the u = 1/2 renormalon: 



5Z V . D , = 5Z P>D , , (28) 

where we have used the fact that the contribution from the upper limit of the integration 
vanishes, since the integral is UV convergent in dimensional regularization. We note that 
diagrams involving transverse gluons, higher-order NRQCD vertices or the gauge fields in 
the D 2 operator all contribute with the wrong power count to yield a pole at u = 1/2. 
(In fact, contributions involving transverse gluons and the part of an V • Aj vertex that 
is proportional to I vanish in the Coulomb gauge.) From (|13|), we obtain the u = 1/2 
renormalon ambiguities in the NRQCD matrix elements: 

A(K. D ,) t = -K±^SL(K.) t , (29a) 

A( Kd ,>, = -*3^«, • (»b) 

Next we identify the u = 1/2 ultraviolet renormalons in the matrix elements (/C)y and 
(fC) p associated with the mixing of these matrix elements into themselves. 
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First consider the heavy-quark and heavy-antiquark self-energy diagrams. It can be seen, 
by counting powers of momentum, that the diagrams involving a transverse gluon or higher- 
order NRQCD vertices do not yield a pole at u = 1/2. The diagrams involving a Coulomb 
gluon contribute to the u = 1/2 renormalon through the renormalization constant of the 
quark (antiquark) wavef unction. For a heavy quark with external energy E and momentum 
p, the Borel transform of the self-energy diagram is given by 

t{E, p) = 4mC F J E + lo _ ip 1 _ iy/2m + u flljhw ' (30) 

The Borel transform of the renormalization constant of the wavefunction is 



dE 



r dH 1 1 

e= p =o ~ F 1 (2^ (l - I 2 /2m + ief I 2 (-P - ie)« ' (31) 



The antiquark self-energy differs from the quark self-energy only in the sign of the Iq term 
in the quark (antiquark) propagator. Hence, these contributions are equal under a change 
of integration variables Iq — > —Iq. It is useful to symmetrize the integrand under this change 
of variables. Then, the sum of the contributions of the quark and antiquark wavefunction 
renormalizations yields the renormalization coefficient 



5Z W { = —AitiCp 




dim dH 



(27T) 



1 

+ 



{Iq - I 2 /(2m) + ie) 2 ' (-/„ - I 2 /(2m) + ief 
l 2 (-P-ie)« ' (32) 



If we carry out the l contribution by contour integration, the individual cut contributions 
yield poles at u = 1/2. However, it is easy to see that the sum of the cut contributions from 
the quark and antiquark terms does not have a pole at u = 1/2. Computing the double-pole 
contribution and integrating over I, we obtain 



^'~3^ {?>) < 33) 

Next we consider the vertex correction involving a Coulomb gluon with A vertices. From 
the right side of Eq. we find that this vertex correction yields the renormalization 

coefficients 




dim dH 1 1 



(2tt) 4 Zo - I 2 /(2m) + it -Iq - I 2 /(2m) + ie 

V(-P- fc )- - (34) 



If we carry out the Iq integration by the contour method, we see that that cut contribution 
does not contain a pole at u = 1/2. However, the contribution of the quark (or antiquark) 
pole is 

167rm ( H 2 \ U f dim d 3 l 1 1 
Mv+ = 5Zp„~ — \^j j x 7^ ?[z2 _j4 /(4m2)K (35) 
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Expanding the last denominator in a power series in l 2 /m 2 , we find that the leading term 
contains a pole at u — —1/2, but not at u — 1/2. The second term in the power-series 
expansion of the last denominator in Eq. fl35|) does yield a pole at u — 1/2: 



R7 R7 AU f^Y f Aim dH 1 

^v >a = SZ P , a ^-—l-j j^j*-^. (36) 



Carrying out the integration over I, we obtain 



u 



i2r ' = t2 "—*L (ft Y^Tu (37) 

There is also a vertex correction involving a Coulomb gluon with Aq vertices and a 
D 4 / (8m 3 ) correction in either the quark or antiquark propagator. The resulting renormal- 
ization coefficients are 

/ fi 2 V r di ™ dH Z 4 1 1 

5z V t = 6z Pfi = 4<iriC F \^) J x j^y 4 (/q - f/{2m) - ie)2 _ /q _ /2/(2m) - ie 

(38) 

where we have doubled the contribution with a D / (8m 3 ) insertion in the quark propagator 
in order to account for the equal contribution with a D 4 /(8m 3 ) insertion in the antiquark 
propagator. Integrating over Iq by closing the contour in the upper half-plane, we obtain 

4:71 ( /i 2 \ " /"dim -y 

5Z v , b = 5Z P ^ — ^-j J x j^^, (39) 

where, as usual, we have dropped the cut contribution, which does not yield a pole at 
u — 1/2. Carrying out the I integration, we identify the u = ~ renormalon: 

^ t = ^_JLg)"_J_ (40 ) 

Next we consider a vertex correction involving a Coulomb gluon with one A vertex and 
one D • E vertex. It yields the renormalization coefficients 

x* xv a n /W ^ m dH l " 1 1 

0Zy c = 0Zp c = ATTtCp I — 



e c J Jx (2tt) 4 4m 2 l - I 2 /(2m) + ie -l - I 2 /(2m) + %e 

X F(^ ' (41) 

where we have taken into account the two diagrams obtained by interchanging the Aq vertex 
and the D • E vertex. This expression differs from ( pB"l) by an overall factor 1/4. Thus from 
(f28|), we obtain the contribution to the u = 1/2 renormalon: 

= ~ sis (42) 
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Finally, we consider the vertex correction involving a transverse gluon with two <x ■ B 
vertices. These vertices violate the spin symmetry, and, so, the vertex correction depends on 
the spin structure of the original matrix element. The resulting renormalization coefficients 
are 

szvAv- AmC /W f dim dH ( ffX| )« 1 - 1 



5Z P)d J \e c J Jx (2vr) 4 2m l - I 2 /(2m) + ie -l - I 2 /(2m) + i 

(<r x l)j Sij - klj/l 2 



2m (-l 2 -ie) 1+u 



(43) 



where 



{'. < 44 > 



Integrating over l by closing the contour in the upper half-plane and simplifying the inte- 
grand by averaging over the angles of I, we obtain 

^} S ~9^^^J A (^FF^' (45) 

where we have discarded the cut contribution, which does not contain a pole at u = 1/2. 
Carrying out the I integration, we can identify the u — ^ renormalon: 



SZ Vtd \ 1 ( H 2 Y 1 r +4/9 



5Z Pd ) urn \e c I 1 — 2u 1—4/3 



(46) 



It may be seen, from power-counting arguments, that the diagrams we have discussed so far 
are the only ones that contribute a pole at u = 1/2. Adding the contributions in Eqs. (j33|), 
(0), (0), (@2|), and (g§) we obtain 



% A_ f/L) 1 J7/9 
SZ P I 7rm \ e c / 1 — 2u 



(47) 



Therefore, from Eq. (|13[), we see that the ambiguities in the operators and JC V are 

A<£*> = -^o^%r^) ' ( 48a ) 
9p e^' z m 

*(Kv) = K J Q ^^n) ■ (48b) 



The results of our calculations of the u — 1/2 renormalon in the operator matrix elements 
for S-wave charmonium decay are summarized in Table ||. 

We note that in NRQCD, in contrast with HQET and with the light-cone expansion, the 
renormalon ambiguities scale differently from the matrix elements themselves. For example, 
the u— 1/2 renormalons computed in this section have ambiguities that scale as Aqcd/^ 
relative to the matrix elements of lowest order in v, but the matrix elements themselves 
scale relatively as v or v 2 . The differences in scaling behavior between the renormalons and 
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the matrix elements arise because the dynamics of heavy-quarkonium system depends on 
the scales mv and mv 2 , rather than the scale Aqcd- 

In Ref. PJ, Braaten and Chen studied the large-order asymptotic behavior of the per- 
turbation series for short- distance coefficients in the NRQCD factorization formulas for the 
decays J/tp —>■ e + e~ and i] c —>■ 77. They calculated the Borel transforms of the short- distance 
coefficients of the leading (in v 2 ) and first subleading matrix elements. By requiring that 
the renormalon ambiguities in the short-distance coefficients be cancelled by the renormalon 
ambiguities in the NRQCD matrix elements, they deduced the u — 1/2 renormalon ambi- 
guities in the matrix elements.^ Comparing the results obtained in Ref. |TJ with the results 
of our direct calculation of the matrix-element ambiguities in Eqs. (|29| ) and P8|), we find 
agreement. Thus, the general principle of the cancellation of renormalon ambiguities be- 
tween the short-distance coefficients and the operator matrix elements is supported by this 
specific example in context of the NRQCD factorization formalism. 



VII. RENORMALONS IN THE GREMM-KAPUSTIN RELATION AND THE 

QUARK MASS 



In this section, we discuss the renormalon content of the matrix elements that appear in 
the Gremm-Kapustin relation WM and the renormalon content of various definitions of the 



heavy-quark mass. 

The Gremm-Kapustin relation is an equality between NRQCD matrix elements that 
relates the heavy-quark kinetic-energy operator to the energy of the quarkonium meson. It is 
correct at leading order in v. One form of the Gremm-Kapustin relation follows immediately 
from the leading-order equations of motion of the NRQCD Lagrangian. In the case of the 
J/ip state, for example, we have 

(0|x f DV^|JM = -m(0|iSb(xW)IW)- (49) 

There are analogous relations for states with different spin and orbital- angular-momentum 
quantum numbers. 

The u — 1/2 renormalon in the matrix element on the left side of Eq. (]49|) is proportional 
to the matrix element of the lowest order operator x The constant of proportionality 
is given by m 2 times Eq. (|28|). The right side of Eq. fl49|) is proportional to the NRQCD 
energy of the quark-antiquark state E$: 

(0|id o (xW)|JM = ^(olxWl J/f)- (50) 

The mixing of id^x^ "^) m t° t ne lower-order operator x^ "^ is obtained by evaluating the 
QCD corrections to the matrix element (0|i<%(x^cn/>)|cc), with the quark and antiquark 



The matrix elements of Braaten and Chen differ from the ones that we consider here by a factor 
\/2Mh, where Mjj is the quarkonium mass. However, Mh is ambiguity free, so this change in 
normalization has no effect on the sizes of the ambiguities relative to the lowest-order-in-t> matrix 
elements. 
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taken on shell and at zero external momentum. Beyond tree level, the mass-shell energy at 
zero momentum is shifted from zero to 



E = E(E,p)\ E=p=0 , (51) 

where T,(E, p) is the quark (or antiquark) self-energy (see Eq. ([30])). This shift in the mass- 
shell position yields a mixing of ido(x^ (T)ip into X^ "^-, with a constant of proportionality 
2Eq. (The factor of two comes from the sum of the equal quark and antiquark energy shifts.) 
Using Eq. ([[4]), we can compute the renormalization coefficient for the mixing of 2Eq into 
unity in passing from the cutoff regulator to dimensional regularization. The contribution 
to the Borel transform of that coefficient that comes from the graph containing a Coulomb 
gluon with two A vertices is given by 



dim rfH 

SZ V:do = 5Z Elp = 25Z Eo = 4:iriC F I 



1 1 

+ 



a (2tt) 4 [l Q -r/2m + ie -l -l /2m + ie 

, (52) 

I 2 H 2 - ie)« 1 ; 

where we have taken into account both the quark and antiquark contributions, which has 
the effect of symmetrizing the integrand under Iq — > —Iq. We have made explicit the fact 
that, to leading order in v, SZv,a = $Z Ei> , since, in Eq. (|50"D, any regulator dependence 
in {Q\x^ cr%l)\J / ip) is suppressed by the factor Ey, which is of order v 2 . The contribution in 
Eq. ( |5"2"D is equal to — m times the contribution in Eq. ( p6|) and yields a pole at u = 1/2: 




SZ V}do = 5Z E ^ = 25Z Eo — ( — ) — . (53) 



1 

37r V, e c ' y 1 — 2u 

Hence, the left and right sides of Eq. ([49]) have equal u = 1/2-renormalon content. 

Now, the NRQCD energy is related to the physical mass of the charmonium state as 

M lP = E,p- 2E + 2m pole , (54) 

where m po i e is the heavy-quark pole mass. In a theory without confinement, such as Quantum 
Electrodynamics, both m po i e and E have nonperturbative definitions in terms of the poles 
in the heavy-quark propagators in the full theory and the effective theory, respectively. In 
a confining theory, such definitions are problematic; only the linear combination m po i — Eq 
has a nonperturbative definition in terms of quarkonium matrix elements. Therefore, we 
absorb Eq into redefinition of m po i e : 

m polc = m pole - Eq, (55) 

which implies that 

2m pole = — E^p. (56) 



Using this redefinition, we can rewrite the Gremm-Kapustin relation ([49]) in the more con- 
ventional form 

(Olx+DVVIJM = -m(M^ - 2m pole )(0|xW| J/^>. (57) 
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In dimensional regularization, Eq vanishes and mp de is equal to m po i e . In the literature, 
m po i e and m pole are used interchangeably. However, it is important to bear in mind the differ- 
ences between m po i e and ^ po i e -0 In a confining theory, m pole is defined only in perturbation 
theory, whereas m' pole is defined, through Eq. (|56|) and Eq. ([50"D, in terms of operator matrix 
elements. In a cutoff scheme, m pole exhibits a dependence, through on the UV cutoff, 
diverging as its first power. (In a dimensional-regularization scheme, this power dependence 
of m' polc on the cutoff disappears.) m' pole also depends on the effective theory in which it is 
defined. In contrast, m po i e is independent of the cutoff and the effective theory. This is true, 
by definition, in a nonconfining theory, where m po i e is the physical mass, and it is enforced 
in perturbation theory in a confining theory. 

In Eq. (|54|) , U V renormalons are absent in Eq and E^ in cutoff regularization and cancel, 
by Eq. ((53|), in dimensional regularization. Since Eq vanishes in dimensional regularization, 
the IR renormalons in Eq must be equal in magnitude and opposite in sign to the UV 
renormalons.Q E^, being an operator matrix element, has no IR renormalons. Therefore, 
we conclude that m po i e has IR renormalon ambiguities, and that near u = 1/2 

4 f^Y 1 

mpole ~ 3^ ) T^Tu- (58) 

Of course, we expect m po i e to contain IR renormalon ambiguities, since the counterterm for 
on-shell mass renormalization involves a loop integration down to zero momentum. It is im- 
portant to recognize that m po i e is neither a short- distance coefficient nor an operator matrix 
element. Rather, m po i e is a quantity that is defined perturbatively in terms of the pole in 
the heavy-quark propagator. Consequently, m pole is factorization-scheme independent, and 
there are renormalon ambiguities in m po i e , even if one employs a cutoff factorization scheme. 
Of course, one could absorb renormalon ambiguities in m po i e into cutoff matrix elements 
of operators in NRQCD. However, the residual renormalon-free short-distance coefficient 
would not be equal to the position of the pole in the heavy-quark propagator. 

We remark that, in a dimensional-regularization factorization scheme, m pole = m pole , 
and, so, 

T^Tu- (59) 

As with m po ie, m-jjs * s neither a short-distance coefficient nor an operator matrix element. 
m I75 i s related to the bare mass m of full QCD as 

m MS = m o + Sm m . (60) 



9 These remarks also are relevant in HQET. A, the heavy-light meson energy in HQET, is analogous 
to E,p in NRQCD. In HQET m pole = M meson — A, where M meson is the physical meson mass. 

10 This is an example of the replacement of IR renormalons with UV renormalons through the 
manipulation of a quantity that vanishes in dimensional regularization. 
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mo is a parameter of full QCD and is independent of any ambiguities associated with factor- 
ization; Sm-jYg consists of poles in e and associated constants, which contain no ambiguities. 
Hence, in contrast with m po i e , m-jjg contains no ambiguities. 

Now we wish to demonstrate that our results for the ambiguities in Ey, Eq and m po i e are 
consistent with the fact that m-jjg has no ambiguities. First we re-express m po i e in terms of 
the heavy-quark MS mass mjjg\ 

Am = m pole - mj^ = S QCD (p)| p . 7=mpolc - Smj^, (61) 

where E^ CD (p) is the heavy-quark self-energy correction in full QCD and (5™-^ is the usual 
MS mass counterterm [Sm-^g = (3CVa s m/47r)(l/e — 7# + m47r) in order a s ]. Thus, we can 
rewrite Eq. ([54]) as 

E i , = My- 2m TTE + (2E - 2 Am). (62) 

Since the expressions for the full QCD heavy-quark energy shift £ QCD (p)| p . 7=mpole and the 
NRQCD energy shift Eq [Eq. fl5T|)l have, by construction, identical behavior at small loop 
momentum, their difference Am — Eq has no support at small loop momentum. Therefore, 
Am — Eq contains no IR renormalon ambiguities. Ey, being an operator matrix element, 
has no IR renormalon ambiguities. Furthermore, Am is free of UV renormalon ambiguities 
because the counterterm 5mj^ removes the poles at u = 0, which correspond to the UV 
renormalons for a logarithmically divergent quantity such as E^ CD . We also see, from 
Eq. (^), that the quantity Ey — 2Eq is free of UV renormalon ambiguities. Thus, in 
Eq. fl62|), mjfg must be free of both UV and IR renormalon ambiguities. 

We note that the absence of renormalon ambiguities in m-jjg implies that Eq. (|58|) is in 
agreement with the standard result |7| for the ambiguity of m po i e relative to rrij^. 

VIII. SUMMARY AND DISCUSSION 

Factorization formalisms for QCD are of great computational utility because they allow 
one to separate short-distance, perturbative physics from long-distance, nonperturbative 
effects. However, the perturbation series for the short- distance coefficients in QCD are, 
at best, asymptotic. Therefore, perturbative calculations of QCD processes are ultimately 
limited in precision. 

One reason for the nonconvergence of the perturbation series is that, even in the case of 
"IR safe" quantities, loop integrations typically extend down to zero momentum. In that 
region, the running coupling becomes large and perturbation theory fails. In the bubble- 
chain model (described in Sec. 0), the coefficients in the perturbation series grow as the 
factorial of the order in a s , leading to renormalon singularities in the Borel transform and 
to ambiguities in the Borel summation of the series. 

By employing a cutoff factorization scheme, in which loop momenta never become small, 
one can eliminate renormalon ambiguities from the perturbation series for the short- distance 
coefficients. We wish to emphasize that the renormalon ambiguities are emblematic of other 
nonperturbative effects that could arise in the low-momentum, long-distance region and 
which are not, as yet, well- understood. A cutoff factorization scheme excludes this nonper- 
turbative region from the coefficient functions. In a dimensional-regularization factorization 
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scheme, one integrates all IR-finite terms in the integrand down to zero momentum. Hence, 
dimensionally-regulated short- distance coefficients contain renormalon ambiguities. 

Physical observables are free of ambiguities, in the sense described in Section |. This 
implies that ambiguities in the factorized expression for a physical observable must cancel 
between the short-distance coefficients and the operator matrix elements (provided that the 
factorization procedure and the underlying effective field theory are sufficiently accurate as 
an expansion in inverse powers of the large scale in the process). That is, the ambiguities 
are merely artifacts of the factorization procedure. 

The presence of ambiguities in the dimensionally-regulated short- distance coefficients 
implies that at least some of the dimensionally-regulated matrix elements contain ambigu- 
ities. Since operator matrix elements are determined completely, in principal, in terms of 
short-distance coefficients and physical observables, the absence of ambiguities in the cutoff 
short-distance coefficients implies that the cutoff matrix elements are free of ambiguities. 
Lattice regularization is one example of a cutoff regulator and, indeed, the QCD operator 
matrix elements are defined unambiguously in lattice regularization. 

The cancellation of ambiguities in a physical observable holds, by construction, if one 
infers the matrix elements through a comparison of theoretical expressions with observables. 
If one determines the matrix elements through a calculation in an underlying effective field 
theory, as we have done in this paper, then the cancellation of ambiguities requires that the 
effective theory reproduce the low-momentum behavior of the original theory. 

We stress that, in computing QCD corrections to the dimensionally-regulated matrix 
elements, it is essential to adopt a convention that is consistent with the standard compu- 
tations of the short- distance coefficients, in which IR-finite expressions are integrated down 
to zero loop momentum. This convention, which is described in detail in Sec. IV B , in- 
volves writing the loop corrections as linear combinations of operators in the underlying 
effective theory and then expanding the loop integrals in powers of the loop momentum 
divided by the large scale in the factorization formalism. The net effect is to remove from 
the matrix elements the power UV divergences that correspond to the IR-finite terms in the 
short- distance coefficients. 

In this paper, we have presented a method for computing the renormalon ambiguities 
in dimensionally-regulated matrix elements. The method exploits the fact that cutoff ma- 
trix elements are unambiguous, which implies that one can compute the ambiguities in 
dimensionally-regulated matrix elements by computing the differences between the cutoff 
and dimensionally-regulated matrix elements. These differences are short- distance quanti- 
ties and, hence, can be computed in perturbation theory. 

We have used our method to compute the u — 1/2 renormalon ambiguities in the matrix 
elements that appear in the NRQCD factorization expressions for the annihilation decays of 
S'-wave heavy quarkonium at the leading and first sub leading orders in v. On comparison 
with the results of Braaten and Chen for the renormalon ambiguities in the corresponding 
short-distance coefficients, we find that the ambiguities cancel in the physical decay rates, as 
expected. The Gremm-Kapustin relation between operator matrix elements is also consistent 
with the ambiguities that we find. Our result for the ambiguity in m po i e is consistent with 
the ambiguity in the expression for m po i e in terms of mjjg given in Ref . || . 

In analyzing renormalon ambiguities, it is important to bear in mind that m po i e and 
m-jjg are neither short-distance coefficients nor operator matrix elements, m-jjg is equal to 
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the sum of the bare mass mo of full QCD and the MS mass counterterm, both of which 
are free of ambiguities and independent of the factorization scheme. m po i e is a quantity 
that is defined in perturbation theory by the on-shell renormalization condition for a heavy 
quark. Because the on-shell renormalization counterterms involve integrations down to zero 
momentum, m pole is ambiguous. m po \ e is also independent of the factorization scheme. 

The accounting of renormalon ambiguities that we have described in this paper can be 
obscured if one writes operator matrix elements in terms of m po i e . For example, the matrix 
element for the energy of the quarkonium state in NRQCD can be written [Eq. ([54])] in 
terms of the physical quarkonium mass M^, m po i e , and the NRQCD energy shift Eq. Since 
Eq vanishes in dimensional regularization, one can use this relation to trade ambiguities in 
the NRQCD quarkonium-energy matrix element for ambiguities in m po i e .Q The analogue of 
this procedure for HQET is frequently employed in discussions of heavy-light mesons. As a 
further step, one can remove the ambiguity in m po i e and the corresponding ambiguity in the 
short- distance coefficients by writing m po i e in terms of m-^jg. 

The renormalization procedure for full QCD can introduce additional IR renormalons 
into the short-distance coefficients if the renormalization counterterms themselves contain 
IR renormalons. This occurs whenever the counterterms involve loop integrations down 
to zero momentum, as is the case in on-shell renormalization. The ambiguities associated 
with these additional renormalons are not cancelled by ambiguities in the operator matrix 
elements, but, rather, by ambiguities in the parameters of full QCD. For example, because 
m po i e is ambiguous, one introduces additional ambiguities into the short-distance coefficients 
if one expresses them in terms of m po i e , rather than, say, nijjg. The additional ambiguities 
in the short- distance coefficients are cancelled by the ambiguities in m po i e . 

In practice, one determines dimensionally-regulated matrix elements either by compar- 
ing perturbative expressions for physical observables with experiment or by computing the 
perturbative relations between dimensionally-regulated matrix elements and some reference 
cutoff matrix elements, such as lattice matrix elements. In both methods, the ambiguities in 
the perturbation series make the determinations of the matrix elements ambiguous. Hence, 
in order to describe a dimensionally-regulated matrix element, one must give not only its 
value, but also the order in ot s in which it was determined. At low orders in a s , this last spec- 
ification may not be important numerically, but at high orders, when the factorial growth 
of the perturbation series dominates, it is essential. 

Finally, let us discuss the relative strengths and weaknesses of the dimensional- 
regularization and cutoff factorization schemes. Dimensionally-regulated short- distance 
coefficients can be computed relatively efficiently, but the computation of, for example, 
lattice-regulated short-distance coefficients, is much more complicated. The dimensionally- 
regulated short- distance coefficients and matrix elements are ambiguous, whereas the cutoff 
short-distance coefficients and matrix elements are not. Of course, the presence of renor- 
malon ambiguities in the dimensional-regularization scheme is not a fatal flaw: In physical 
observables, the factorial growth in the dimensionally-regulated short- distance coefficients 
cancels against a similar growth in the dimensionally-regulated matrix elements. However, in 



This relation, with Eq set to zero, is often taken as a definition of the pole mass. However, one 
can neglect Eq only in a dimensional-regularization factorization scheme. 
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order to achieve this cancellation, one must work to the same order in a s in the short-distance 
coefficients and the matrix elements. In addition, it is essential to work to sufficient accuracy 
in the small scale of the underlying effective field theory JT0[ . In practice, one cannot achieve 



a complete cancellation, but one can systematically reduce the renormalon ambiguities by 
introducing more operator matrix elements and/or including more terms in the effective ac- 
tion. In general, cutoff matrix elements are very sensitive to the cutoff (factorization scale), 
typically scaling as a power of the cutoff. This sensitivity to the cutoff is cancelled by a 
similar dependence in the short-distance coefficients, but the cancellation is imperfect at any 
finite order in perturbation theory. A related difficulty is that the cutoff matrix elements 
exhibit their nominal sizes only if the cutoff is chosen to be of the order of the small scale in 
the calculation. In contrast, dimensionally-regulated (MS) matrix elements have no power 
dependence on the cutoff. Nevertheless, at high orders in ac s , they must deviate from their 
nominal sizes because of the factorial growth of the determining perturbation series. 

Of course, in the end, the physical observables are independent of the choice of factoriza- 
tion scheme, provided that one works to sufficient accuracy in the small scale of the effective 
theory and to sufficient and consistent accuracy in a s . However, it is important to bear 
in mind that the properties of the dimensionally-regulated and cutoff quantities are rather 
different and that the differences may be significant in practical calculations. 
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TABLES 



Diagrammatic Element 

quark propagator 

Aq vertex 

V • Aj vertex 

Aj • Aj seagull vertex 

cr • B spatial-gluon vertex 

D • E temporal-gluon vertex 

D • E spatial-gluon vertex 

D x E • a temporal-gluon vertex 

D x E • <x spatial-gluon vertex 

D x E • er spatial-temporal seagull vertex 

D x E • a spatial-spatial seagull vertex 

D 4 /(8m 3 ) quark-propagator correction 

D 4 /(8m 3 ) spatial-gluon vertex 

D 4 /(8m 3 ) spatial-gluon seagull vertex 

D 4 /(8m 3 ) 3-spatial-gluon vertex 

D 4 /(8m 3 ) 4-spatial-gluon vertex 



Feynman Rule 
i/[±Po - p 2 /(2m) + ie] 

ig(Pi +P'i)tal{^rn) 
-ig 2 (5ij t b t a + perm)/ (2m) 
ge ijk hj(Tk t a /(2m) 
±igl 2 1 t a /(8m 2 ) 
^igluho t a /(8m 2 ) 
±gtijkPiPjVk t a /{4:m 2 ) 
Tgtijkip + P')j0kho t a /(8m 2 ) 
±g 2 (eijkhjCTk tbt a + perm)/(4m 2 ) 
±g 2 (eijk<7kho tbt a + perm)/ (4m 2 ) 
ip 4 /(8m 3 ) 

-ig(p 2 + P' 2 )(p + P')i t a /(8m 3 ) 

ig 2 {[(2p> - l 2 ) j (2p + h)i + (p 2 + p' 2 )5 i3 } t b ta 

+ perm} /(8m 3 ) 
-ig 3 {[(2p' - l 3 ) k 6ij + (2p + h)i6 k j] t c t b t a 

+ perm} /(8m 3 ) 
ig 4 (Sij6km t d t c t h ta + perm)/ (8m 3 ) 



TABLE I. Feynman rules for quark and antiquark propagators and interaction vertices in 
NRQCD. The upper (lower) signs correspond to quarks (antiquarks). We have retained only the 
"Abelian" terms, i.e., those that contain no color-matrix commutators. Energy and momentum 
are conserved at the vertices, and we observe the convention in which the antiquark diagrammatic 
energy and momentum are the negatives of the physical energy and momentum, p is the incom- 
ing (diagrammatic) quark or antiquark momentum, p' is the outgoing (diagrammatic) quark or 
antiquark momentum. The gluons have incoming momenta, spatial polarization indices, and color 
indices a), (l2,j,b), (l^,k,c), and (l^,m,d), respectively. T denotes an SU(3) color matrix in 
the fundamental representation, normalized to tvt a tb = (l/2)5 a b- "Perm" denotes all additional 
permutations obtained by interchanging the momentum, polarization, and color labels of one gluon 
with those of another gluon. 
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Quantity 


Description 


Spin Triplet 


Spin Singlet 


Correction to {K? D ): 








<* Z V,D 2 /S Z P,D 2 


Aq-Aq vertex correction 


+8/3 


+8/3 



Corrections to (K): 

a 

SZ Vib /5Z Ptb 
5Z v ,c/$Zp )C 
$Zv,d/5Zp jd 
5Z V /SZp 



Quark wave- function renormalization +2/3 +2/3 

Aq-Aq vertex correction — 1/3 — 1/3 

Aq-Aq vertex correction with D 4 prop, insertion —2/3 —2/3 

A -Y) ■ E vertex correction +2/3 +2/3 

<r • B-er • B vertex correction +4/9 —4/3 

Total correction to (K.) +7/9 -1 



TABLE II. The u = 1/2-renormalon contributions in the mixing of the spin-triplet matrix 
and spin-singlet matrix elements {K, 2 D )y and (JCfyp into the spin-triplet and spin-singlet matrix 
elements {JC)y and (/C)p, respectively, and the u = 1/2-renormalon contributions in the multi- 
plicative renormalization of the spin-triplet and spin-singlet matrix elements {JC)y and (IC)p. The 
values displayed are the coefficients of (l/7r)(/i 2 /e c *)"(l — 2-u) -1 times the matrix elements. 
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